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Abstract 

Neutrinos with a magnetic moment \i change their helicity when interacting with 
an electromagnetic field. Various aspects of this effect have been described as spin 
precession, spin-flip scattering, and magnetic Cherenkov radiation. These perspec- 
tives are unified in an expression for the vl — * transition rate which involves the 
correlators of the electromagnetic field distribution. Our general formula corrects a 
previous result and generalizes it to the case where the fields cannot be viewed as 
classical and where the momentum transfers need not be small. We evaluate our re- 
sult explicitly for a relativistic QED plasma and determine the depolarization rate to 
leading order in the fine structure constant. Assuming that big-bang nucleosynthesis 
constraints do not allow a right-handed neutrino in equilibrium we derive the limit 
fi < 6.2xl0 -11 hb ° n the neutrino magnetic moment. Bounds on fi from a possible 
large scale magnetic fields are found to be more stringent even for very weak fields. 
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1 Introduction 

A neutrino or other neutral particle with a magnetic moment fi gets depolarized when 
traversing a random distribution of electromagnetic fields as, for example, in a plasma of 
charged particles. This effect could be important in stars or in the early universe where the 
standard weak interactions produce only left-handed neutrinos. More than twenty years 
ago the depolarization effect was considered as a possibility to solve the solar neutrino 
problem, but for plausible values of /i the rate was found to be too small J1J. In supernova 
cores, the left-handed neutrinos are trapped while the helicity-flipped states could escape 
more easily because they scatter only by the assumed magnetic-dipole interaction. The 
depolarization effect would thus disturb the standard picture of supernova theory so that 
the known properties of supernovae and the observed SN 1987A neutrino signal allow one to 
derive certain limits on Dirac magnetic moments |2|, [3|. In the early universe, the radiation 
density and thus expansion rate would increase due to the new thermally excited degree 
of freedom, modifying the predicted primordial light-element abundances. A comparison 
with observations again allows one to derive limits on \x [|], ^ 

In most of these studies the depolarization rate was calculated from the spin-flip scat- 
tering process vl + X — > X + v R where X represents electrons, positrons, or other charged 
particles and the interaction is due to the assumed neutrino magnetic moment. In addition 
there may be electromagnetic modes with a "space-like" dispersion relation which allow for 
the Cherenkov emission v L — > u R + 7 and absorption 7 + v L — > v R and thus contribute to 
the depolarization rate |], || . This is the case in supernova cores where the photon disper- 
sion relation is dominated by the nucleon magnetic moments or in nonrelativistic plasmas, 
but it is not the case in the relativistic e + e~ plasma of the early universe. The rate for 
higher-order scattering processes such as magnetic Compton scattering 7 + v — > v + 7 
is proportional to fi A and is thus ignored in the present discussion which is limited to [i 2 
effects. To order fi 2 right-handed neutrinos can also be produced by pair processes of the 
type 7 vlVr (plasmon decay) or e + e~ — ► vl^r (pair annihilation). In supernova cores or 
the early universe with a large population of left-handed neutrinos the Vl —>■ v R reactions 
always seem to be more important than the pair processes which, however, dominate in 
normal stars which do not have a population of trapped left-handed neutrinos. 

Loeb and Stodolsky || were the first to look at the depolarization effect from a more 
classical perspective. They argued that in a macroscopic magnetic field, left-handed neu- 
trinos simply spin-precess and that even in a microscopic distribution of random fields 
the depolarization rate should be calculable as a suitable ensemble average over the spin- 
precession formula. They found a result which represents the depolarization rate in terms 
of certain correlation functions of the electromagnetic fields. In principle, this correlator 
approach incorporates all electromagnetic effects to order fi 2 which lead to depolarization 
such as spin-flip scattering and the Cherenkov processes. 

These different approaches have co-existed in the literature with no apparent attempt to 
compare them, to understand their relationship, or to check their mutual consistency. One 
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problem is that Loeb and Stodolsky's approach is entirely classical so that it is not obvious 
if and how their formula is limited when applied to a plasma which involves nonclassical 
(quantum) excitations of the electromagnetic field. 

We thus reconsider neutrino depolarization both classically in the spirit of Loeb and 
Stodolsky (Sec. |2|) and from a quantum-kinetic perspective (Sec. [$]). Put another way, we 
derive the depolarization rate in terms of correlation functions of the electromagnetic fields 
which may be quasi-classical or true quantum variables. Even in the classical limit our 
general result involves terms which are absent in Loeb and Stodolsky's formula. 

We also compare with the imaginary part of the neutrino self-energy (Sec. |j) which has 
a simple relation to the depolarization rate. This calculation amounts to a determination 
of the dominant spin-flip scattering rate, but even to lowest order in the fine-structure con- 
stant a the correct screening prescription has to be incorporated by resumming the photon 
propagator. Then it automatically includes all order /i 2 contributions to the production 
rate of right-handed neutrinos, i.e. it includes Ui — > z/r transitions from spin-flip scattering 
or from Cherenkov processes as well as plasmon decay and e + e~ annihilation. 

In Sec. [5] we evaluate the depolarization rate explicitly by using the correlators for a 
relativistic QED plasma and compare the classical and quantum treatments. In Sec. ||] we 
apply this result to the depolarization of neutrinos in the early universe, allowing us to 
derive a limit on the magnetic dipole moment from considerations of big-bang nucleosyn- 
thesis. We also compare this limit to previous bounds obtained from assuming a large-scale 
background magnetic field. In Sec. [7| we summarize our findings. 

2 Classical Trajectory Approach 

2.1 Depolarization Rate 

A nonrelativistic neutral particle with a magnetic moment \i and an intrinsic angular 
momentum (spin) precesses in the presence of an external magnetic field B according 
to P = <7/iB x P. Here, P is the spin polarization vector which, for a spin-| particle, 
parametrizes the spin density matrix in the usual form p = ~(1 + PjCTj). Further, g is the 
gyromagnetic ratio. Loeb and Stodolsky used g — 1, i.e. a classical particle, while we 
always study neutrinos with spin ~ so that g = 2. If the particle moves relativistically, 
the main modification is that in the laboratory frame only the magnetic field transverse to 
the direction of motion contributes, and that a transverse electric field is also important 
because in the particle's rest frame it manifests itself as a magnetic field. Altogether, the 
precession formula for an ultrarelativistic {v = 1) spin-| particle is 

P = 2fi(B ± -pxEjxP , (2.1) 

where the subscript _L refers to the field vectors transverse to the direction of motion and 
p is a unit vector in the direction of the neutrino momentum. 
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The same equation of motion is obtained if one begins directly with the covariant 
Lagrangian which describes the coupling of a Dirac neutrino with magnetic moment /i to 
the electromagnetic field tensor F^ u , 

Ant = -|^^F^* = /i*(S-B-za-E)* , (2.2) 

where \1/ is the neutrino Dirac field. In terms of the usual Dirac matrices we have Ej = 
757°7i and oti = 7°7j. This Lagrangian shows that the magnetic- moment interaction indeed 
only couples neutrinos of opposite chirality. As we are concerned only with ultrarelativistic 
neutrinos this is tantamount to a coupling of opposite helicities. 

In the presence of other forces than those represented by the electromagnetic fields, the 
off-diagonal elements of the density matrix can be damped by collisions which "measure" 
the helicity content of a neutrino state. Standard weak interactions would have this prop- 
erty because only left-handed neutrinos are affected by collisions. This damping effect is 
important when the spin precesses in macroscopic magnetic fields f?|, H , and can be taken 
into account by adding a term — D Pi to the right hand side of the equation of motion for 
Pi, i = 1,2, in Eq. ( |2.1| ). For depolarization in stochastic electromagnetic fields one can 
speak of a coherent spin precession only for a time period which represents the correlation 
time of the fields. The damping effect would be important if the average time between 
weak collisions were shorter than this coherence time scale. An example is a putative 
horizon-scale magnetic field in the early universe. 

We always assume ultrarelativistic (i.e. effectively massless) neutrinos. However, in 
a medium the dispersion relations for left- and right-handed states are different because 
only the left-handed ones feel the weak potential produced by the background medium. We 
may write the refractive energy difference between left- and right-handed states in the form 
u> Te f T = 2fiB Te f r in terms of an effective magnetic field -B re f r which points in the neutrino's 
direction of motion. Therefore, the precession formula P = 2/iB eff x P finally involves an 
effective magnetic field with H c f = — p x Ej_ and Bjj ff = (a; re fr/2/i)p. 

In general, the electric and magnetic fields depend on location and time in arbitrary 
ways. Therefore, a neutrino with a magnetic moment will be deflected. However, the 
deflecting forces on a magnetic dipole are proportional to the field gradients while the 
precession effect depends on the fields directly. Therefore, if the spatial variations are 
small the neutrinos can still be assumed to move on a straight line which can be taken 
to be the z-direction. Moreover, the neutrino can be taken to "see" only the fields at 
a specific location which is assumed to vary with time as z = t because of the assumed 
propagation with the speed of light and because we take z = at t = 0. Therefore, through 
the condition z = t the fields B efr are to be viewed as functions of time alone. We call this 
the "classical trajectory approach" to the problem of neutrino spin depolarization. 

In order to derive the equation of motion of the polarization vector in a stochastic field 
distribution it proves useful to write the precession equation in the form 



P(t) = M(t) P(t) 



(2.3) 
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where the time-dependent matrix M is explicitly given by 



Mi. 



—D 

0> re f r 



E, 



— Wrefr 2/i(_B y — E x ) ' 

-£> -2/i( J B a: + E y ) 
2fi(B x + E y ) o y 



(2.4) 



if the neutrinos are ultrarelativistic and move in the ^-direction. 

In general we must view the matrix M as consisting of one part (M) which is nonzero 
on the average plus stochastic fluctuations around that average. Even if there are no large- 
scale magnetic fields, the refractive effect provides a nonvanishing average contribution. 
In order to derive the depolarization rate we eliminate (M) by going to an "interaction 
picture" with Q(t) = e^ M ^P(t) so that we are left with the equation of motion 

Q(t) = m(t)Q(t) (2.5) 

with 

m(t) = e- (M)t (M - (M))e {M)t . (2.6) 
The formal solution to Eq. ( |2.5| ) is 

°° ft rtn-l 

Q(t) = E / ■ ■ ■ / dt n mfc) ■ ■ ■ m{t n ) Q(0) . (2.7) 
n=0 Jo Jo 



This is the solution we need to average over a statistical ensemble of field configurations. 
We shall assume that the field fluctuations are Gaussian so that the n-point correlation 
functions can be reduced to products of pair correlation functions. Let us observe that in 
QED the high temperature effective action is in fact Gaussian 0. (This is not a general 
feature and not true e.g. for QCD.) We seek the solution only for times t much larger than 
the correlation time of the stochastic fields. In the integral in Eq. ( |2.7| ) all time arguments 
have to occur in close pairs for the integrand to be non- negligible. Since the time arguments 
are ordered it is only the permutation where adjacent matrices are contracted that gives 
any contribution in the leading large-time limit. It can also be checked at the end that this 
asymptotic region is reached long before one decay time in the present case of neutrino 
spin oscillations. We can therefore approximate (Q(t)) by 

°° ft rt n - l 

(Q(*)> ~ E / <tti-~ dtn {m{tx)m{t 2 )) ■■■ (m(t n _i)m(t„)) Q(0) . (2.8) 
^ Jo Jo 



After computing the pair correlation matrix (m(ti)m(t2)) ^ with zero average field strength 
we find that there is no mixing between the third component and the rest so that it is useful 
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to define 
n(ti - 1 2 ) 



-(m(t 1 )m(t 2 )) 33 (2.9) 
(2 f i) 2 e- D ^-^i [ cos[cu Icb (ti - t 2 )](B ± (t 1 ) ■ B ± (t 2 ) + E ± (t x ) ■ E ± (t 2 ) 



+ P 



B(ti) x E(* a ) - E(ti) x B(t 2 ) 



+ Sin^efr^ - t 2 )](E ± (t 1 ) ■ B ± (t 2 ) - Bj_(tl) • E ± (t 2 ) 



p • 



B(ti) x B(t 2 ) -E(ti) x E(t 2 ) 



Of course, in the terms involving cross products of fields we could have equally used the 
transverse field components. 

The depolarization rate is given by the shrinking rate of (Ps{t)). If we carry out all 
the integrals in Eq. ( |2.8| ), keeping only the leading term at large t, we finally find that the 
ensemble-averaged evolution of P 3 is given by 



(P 3 (*)> = exp 



-t I dt'nit') 
h 



P 3 (0) = e" rdc p° 1 *P 3 (0) 



(2.10) 



If u; re f r is much larger than the inverse correlation time of the electromagnetic fields, the 
oscillating cos(u; r ef r i) and sin(u; re f r t) terms under the integral in Eq. (|2.9| ) would suppress 
the depolarization effect. Likewise, a very large helicity-measuring damping coefficient 
D would prevent oscillations and suppress the depolarization. In stars, even supernovae, 
and in the early universe these suppression effects are insignificant because the inverse 
of the electromagnetic correlation length is much larger than both the refractive energy 
difference between those states and the damping rate. Therefore, for all situations which 
are of astrophysical interest we can put u; re f r = D — 0. (In the presence of a large-scale 
background field there is an additional component to Tdepoi from (M) causing a coherent 
spin precession. Of course, this additional term would depend on both u re f r and D; see 
Sec. |Q| .) The depolarization rate from stochastic fields is then found to be 

poo . 

r dep oi = (2/i) 2 ^ dt(B ± (t)-B ± (0) + E ± (t)-E ± (0) 

+ p- [B(t) x E(0) -E(t) x B(O)l) . (2.11) 



This result agrees with that of Loeb and Stodolsky || except for the cross term. They 
agree that it should be there and stress that it can be derived rather easily by beginning 
with the spin-precession formula in the neutrino's rest frame where only the B^ correlator 
appears, and express it by the Lorentz transformed laboratory fields.Q 

^^L. Stodolsky, private communication. 
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Two remarks about the interpretation of Eq. (|2.11 ) are in order. First, a term like B(£) 
really means B(£, r) with r = pt where p is the neutrino velocity vector. Therefore, even 
in an isotropic medium a vectorial quantity like (E(i) x B(0)) need not vanish because it 
depends on the external vector p. Second, the integrand is an even function of t because 
the fields are classical variables so that their ordering is arbitrary, and because of time 
translational invariance of the statistical ensemble in equilibrium. Therefore, the time 
integral may be extended to — oo if a compensating factor 1/2 is introduced. We will 
always use this more symmetric form which is easier to handle in Fourier space, and which 
allows for a direct transition to the quantum case (Sec. ^) where the ordering of the fields 
is important. 



2.2 Isotropic Medium 

In a homogeneous and isotropic plasma or other medium, correlator expressions like the 
ones appearing in Eq. ( |2.11| ) are conveniently calculated in Fourier space. To this end we 
introduce the usual notation^ 

/+oo 
dfe'^^t, k)tyO,-k)) (2.12) 
-oo 



where X\ and Yj each represent a component of the E- or B-fields with Xi(t, k) the spatial 
Fourier transform of Xi(t,r) and so forth. Therefore, the depolarization rate Eq. ( |2.11| ) 
can be expressed as 

r — ¥£>/ S/r*-*—^- <- 3 » 



where the P-dependent "dynamical structure function" is 

S P (K) = (B^ + El + p- (B x E-E x B)) . (2.14) 



It will turn out that fi 2 Sp(K) is to be interpreted as the probability for a neutrino of four 
momentum P to transfer the four momentum K to the medium. The dt integration yields 
27r5(p • k — k ) which allows us to perform the dko integration trivially. Altogether we find 

r d 3 ~k 

r dep oi = 2/i 2 / — - S P (K) , (2.15) 



where K is restricted by the condition ko = p • k. 

The structure of Sp still depends on the field components transverse to the neutrino 
momentum. We stress that contrary to Ref. || the assumed isotropy of the medium does 

2 We use K = (fco,k), P = (j>o,p) etc. for four-vectors and k = |k|, p = |p| etc. for the corresponding 
three- vectors. 
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not imply that (B±)k equals ^{B 2 ) K because (B±)k depends on the external vector p and 
thus is not a scalar. In order to use isotropy properly we observe that B T = B — (p • B)p 
so that 

S P (K) = {B t B j + EiEj) K fa - ptpj) + (BiEj - E^k e m h . (2.16) 

We further note that in an isotropic medium any correlator expression of the form (XjYj)K 
can be proportional only to Sij, kikj, or e^ke because only the vector k is available to 
construct spatial tensor structures. The most general structure is thus found to be 

(X t Y 3 ) K = (X ■ Y) K *S ~ + ((X ■ k)(Y • k)) K 3kikj ~ 6tJ 

+ (XxY) r k^ . (2.17) 
It is then straightforward to show that 

S P (K) = (B 2 + EV 1 + ^ ' + ((B ■ k) 2 + (E ■ k) 2 ) K 1 ~ 3( 2 k ' 

+ (k-p)(BxE-ExB) r k . (2.18) 

This expression simplifies further if we observe that E and B must obey Maxwell's equa- 
tions. Because the magnetic field is divergence free we have B k = 0. Further, ((E-k) 2 )^ = 
(E 2 k 2 - (k x E) 2 )x- Maxwell's equations give us k x E = A; B so that ((E • k) 2 ) K = 
(E 2 )k - (k /k) 2 (B 2 ) K . Next, (B x E) K • k = (B • (E x k))^, allowing us again to apply 
k x E = A; B so that (B x E - E x B) K ■ k = -2(k /k) (B 2 ) K . Altogether we thus find 

S P(K) = (BV (1±C|^! - % p . t - + <ev [i - (P ■ if 

(2.19) 

In an isotropic medium the (. . )k expressions depend only on (k , k) and not on the 
direction k. 

In order to present our final result we note that the depolarization rate studied thus 
far represents the rate by which a fixed ensemble of neutrinos gets depolarized. In the 
early universe or in supernovae, however, a more relevant quantity is the spin-flip rate 
which measures the speed by which the sea of right-handed neutrinos gets populated if 
the number density of left-handed neutrinos is held fixed at its thermal equilibrium value 
because they are replenished by other reactions. Evidently, since the number of right- and 
left-handed neutrinos is = |(1 ± P 3 ) we have that n R /riL = r flip = r dcpol /2 so that 



i _2 
flip — H 



d 3 k 

(2^p 



l-(p.k) 2 l(W+ 1 3( f k)2 (B%) , (2.20) 



where, again, K is constrained by k = p • k. 
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3 Quantum Kinetic Approach 

3.1 Relaxation Rate 

The spin-flip rate derived by the classical-trajectory approach in the previous section is 
valid only in the approximation that the momentum transfer k is always much smaller 
than the neutrino momentum p so that it is justified to view the neutrino as propagating 
on a straight line during a typical field correlation time. Further, the electric and magnetic 
fields were taken to be classical variables, ignoring the quantized nature of the exchange of 
energy and momentum (ko, k) between the neutrino and the medium. 

In a quantum-kinetic approach the relaxation of the neutrino helicity by electromagnetic 
interactions is described by a Boltzmann collision equation of the form 

dtf« = J ^ [WWQ, P)/ q L (l - /*) - W R ^ L (P, Q)f*(l - / q L ) + • • • ] (3.1) 



where f£ ,R are the occupation numbers for left- and right-handed neutrinos of momentum 
p, respectively. Here, Wr-*l(P, Q) is the transition rate for a right-handed neutrino of four 
momentum P to a left-handed one of Q under the influence of the electromagnetic fields of 
the medium. There are also terms representing pair production and annihilation processes 
which we will discuss below. If there are large-scale magnetic fields the Boltzmann equation 
includes an oscillation term analogous to the simultaneous treatment of collisions and flavor 
oscillations in Refs. [10|. 



The neutrino interaction with the electromagnetic fields is of the current- current form. 
It is well known from linear-response theory that in this situation the transition probability 
is given essentially by the dynamical structure function of the medium (in our case the 
electromagnetic fields), i.e. by correlator expressions like the classical ones discussed in the 
previous section. Therefore, all that remains to be done to derive W(P, Q) is to determine 
the required tensorial contraction between the neutrino momenta and the electromagnetic 
field correlators. 

To this end we begin with the interaction Lagrangian Eq. ( j2.2| ) and consider the matrix 
element for the transition of a left-handed neutrino with four momentum P to a right- 
handed one with Q. With the massless neutrino Dirac spinors Up 7 l and uq.r which involve 
the appropriate chirality projections one finds M. = {lUp^H ■ £ — zE ■ ot)uQ t R. Of course, 
because the interaction couples only states of opposite chirality it would have been enough 
to include one chirality projector. Taking the square of this matrix element and carrying 
out the Dirac traces we find 



+ (BiEj — EiBj) K eiji 



2 

Pi + qt 
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+ i(E i E j ) K e ijl ^^ , (3.2) 



where K = P — Q. In the limit of small k we may set p = q, taking us back to Eq. ( |2.16| ) 



of the classical-trajectory approach apart from the new term which involves (E x E)#. It 
is easy to show that in the classical limit where the fields are commuting variables this 
term disappears under the d 3 k phase-space integration so that our present result leads to 
the same classical spin-flip rate. 

In the general case the status of the E x E term is more subtle. It has the opposite 
sign for an R — > L transition. On the other hand it is easy to show that Wl-+r(P, Q) = 
Wr^l{P,Q) if the medium is invariant under parity. Therefore, in such a medium the 
(ExE)jf term must vanish even in the general case. This is not surprising in the sense that 
the electromagnetic properties of a parity invariant medium are characterized by precisely 
two independent "material constants" which can be chosen to be the dielectric permittivity 
and the magnetic permeability. Another pair of equivalent parameters are the longitudinal 
and transverse polarization functions Ht,l{K) which we will use below. Therefore, there are 
only two independent field correlator expressions, e.g. (E 2 ) K and (B 2 )^. They are related 
to Ut,l{K) by the fluctuation-dissipation theorem. A parity- noninvariant medium, on 
the other hand, is characterized by three independent parameters — there are two different 
transverse polarization functions. The third function gives rise to a nonvanishing (E x E) k 
field correlator. 

For the neutrino spin relaxation problem we inevitably have a medium which is not 
parity invariant because initially it involves only left-handed neutrinos. By their assumed 
magnetic moments they are expected to give rise to a nonvanishing (E x E) term, i.e. 
L — > R and R — > L collisions are not expected to occur with the same transition probability. 
However, in terms of the relaxation rate this would be an order /z 4 effect so that to order /i 2 
we may ignore electromagnetic neutrino-neutrino interactions. In the cases of interest to 
us we may thus ignore the (E x E) K term and use W(P, Q) = W L ^ R (P, Q) = W R ^ L (P, Q) 
for the transition rate. 

We proceed by applying Maxwell's equations in the same way as in the previous section, 
never changing the order in which the non-commuting field operators appear. Since we 
never changed this order even in the classical case, we find the same simplifications as 
there. The final contraction of indices then yields 



fi- 2 W(P,Q) = (B )k 



l + (p-k)(q-k) 2fc p-k + q-k fc 2 p-q-3(p-k)(q-k)~ 



k 2 k 2 




This complicated looking expression can be transformed to 



8k 2 p q 



( Ik 2 K 2 
2(EVl(B\l4 + 



k 2 (Po + qo) 2 , 



(3.3) 



(3.4) 
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Note that q = p — k and p = p and that — K 2 = — (P — Q) 2 = 2PQ = 2(p q — p • q) is 
a positive number. 

We now return to the Boltzmann collision equation ( |3.1| ) for a parity invariant medium 
where we need only one function W(P,Q). The transition rate for the pair production 
and annihilation processes is given by the same function with "crossed" four momenta, i.e. 
P — > —P or Q — > —Q. The complete collision equation is then 

dtf? = J^[w(Q,P)rt(l-f*)-W(P,Q)f«(l-ft;) 

+ W(-Q,P)(l-rt)(l-fK)-W(Q,-P)fZf*} . (3.5) 

Here, the first term in the collision integral represents gain by L — > R spin-flip scattering, 
the second loss by R — > L scattering, the third gain by pair production, and the fourth 
loss by pair annihilation. 

In the astrophysical applications that we are interested in (supernovae, early universe) 
the left-handed neutrinos are and remain in thermal equilibrium so that we may replace f£ 
by the Fermi-Dirac distribution f£ at the appropriate temperature and chemical potential. 
The collision equation is then a linear differential equation of the form 

dig = r g ain( P )(l " /*) " r ioss ( P )/ p R (3.6) 



with 



r g ain(p) 



rioss(p) 



d 3 q 

(2^)3 

d 3 q 

(2tt) 3 



W(Q,P)fF + W(-Q,P)(l-fF 



W(P,Q)(l-e) + W(Q,-P)f> 



Further, in equilibrium we have d t f^ = and = f£. Inserting this in Eq. (|3] 



that r gain = r to t/p with r to t = r gain + ri oss so that we may write 



(3.7) 
reveals 



W P * - ft 



-rtotCpJC/?-/, 



{3.1 



Therefore, r tot (p) is the appropriate rate that measures the exponential approach of the 
p mode to helicity equilibrium. 



3.2 Correlation Functions 

In order to evaluate the relaxation rate we need to know the electric and magnetic field 
fluctuations in a given medium. By virtue of the fluctuation-dissipation theorem they are 
related to the imaginary part of the medium's dielectric response functions. One way of 
expressing these quantities is in terms of the longitudinal and transverse photon spectral 
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functions At,l which are the coefficients in the decomposition A^ v = —V^Ar — Q^uAl 
in the Landau gauge. The photon spectral function A^(K) is related to the retarded and 
advanced Green's functions through A^K) = [iD^ u (ko + «e, k) — iD^ u (ko — ie, k)]/27r. 
The transverse and longitudinal projection operators have the nonvanishing components 
Vij(K) = -5 i:j + (hk^/k 2 and Q^{K) = — (K 2 k 2 )~ l (k 2 , — A; k) At (A; 2 , — /c k)„. In terms of 
the photon polarization functions we have 

1 Im Tlx i 



AtAK) = - 



71 



|fr 2 -Ren TiL | 2 + |imn TiL | 2 



(3.9) 



The analytic continuation in the imaginary part is the retarded ImlLr^/co + ^e). Note that 
the photon polarization tensor is given by = V^Ht + Q^v^l- 



In terms of the spectral functions the field fluctuations are found to be [ITT 

2tt 



1 _ e -/3fco 



(E 2 ) 



2tt 



K 



2k At(K) 



2k 2 A T (K) + K 2 A L {K) 



(3.10) 



1 _ e -/3fco 

where /3 — 1/T is the inverse temperature. Note that for a positive ko (energy given to 
the medium) the thermal factor is identical with 1 + (e /3fe ° — l) -1 , i.e. it is understood as a 
Bose stimulation factor for exciting a quantum (ko, k) of the medium. Conversely, if k < 
(energy lost by the medium) this factor is — (e^ — l) -1 . Apart from the minus sign it 
is the occupation number for such an excitation. The functions At,l(K) are odd in ko so 
that the negative sign of the thermal factor is automatically absorbed. Put another way, 
the correlators obey detailed-balance conditions of the form (B 2 )_^ = (B 2 )ft-e~^ fe °. 



4 Imaginary Part of the Neutrino Self-Energy 



An alternative way of calculating the rate of populating right-handed neutrinos is through 
the imaginary part of the neutrino self-energy [|TI|. The relaxation rate L tot introduced in 
the previous section is directly related to the imaginary part of the neutrino self-energy 
through 

2lmu P)R 'L(po + ie)u PjR 



r t ot(p) 



Up,RU PtR 



(4.1) 



where Up y p is the Dirac spinor of a right-handed neutrino with four momentum P. 

The self-energy to one-loop order is a bubble diagram with a left-handed neutrino 
and a photon propagator connected via the interaction vertex fiK a a afJ- . Its time ordered 



imaginary part can be shown to be JL3 

Afao) 



ImS(P) 



sin 2<p p 



xK n a 



d A K . , . .1 1 , „ 

j^^e(p + ko)e(ko)-sm2(pp +K -smh29 K 

>+ /C)i(l - l5 )K^(2n) 2 S((P + K) 2 )A,„(K) 



(4.2) 
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where e(x) = ±1 depending on the sign of x, the photon spectral function A^ u was defined 
in Sec. |372| above, and 



- sin 2(j> K 



e P\ko\/2 
of\k \ _|_ I 



- sinh 29 k 



oP\ko\/2 



e /3|fco| 



(4.3) 



The neutrino inside the loop is necessarily left handed if the external one is right handed 
since the interaction flips chirality, so we do not need to include the |(1 + 75) part of the 
propagator to compute the right-handed self-energy. The fermion distribution function in 
Eq. ( |4.3|) should therefore describe a fully populated equilibrium ensemble for which the 
chemical potential is taken to vanish. 

The results of Refs. |12| and |13] differ by the overall factor e(po) which is related to 
the retarded [E| or time ordered |13| prescription for ImE. The physical sign is however 
obvious since r tot must be positive so we can concentrate on po > 0. Explicitly we find for 
an on-shell vr with momentum p 



r to t(p) = 7t\ dkk dk e(-K 2 (K 2 + Apk + Ap 2 
Ztt Jo J -00 \ 



0) 



o p\p +k \ 



+ — — + e(p + k )8(-k ) 



,0\ko\ _ 1 



e(k )e(-p - k ) 



x- 



K 



1 H — + -7—9 M*Q ~ 1 + 7T MK) 



p 4p 



2p 



e(k 



(4.4) 



where the medium was taken to be parity invariant so that there is only one transverse 
polarization function. There are contributions to r tot both from creation and annihilation 
of vr through several processes. Depending on the signs of k^ and p+ko these processes can 
be divided into pair creation/annihilation, Cherenkov radiation and scattering off charged 



particles through virtual photon exchange [JT2] . Calculating the imaginary part of the self- 
energy is thus a convenient way of obtaining all processes that are allowed at the one-loop 
level with a interacting photon correlation function. 

As expected, this relaxation rate is identical to r gain + T\ oss from Eq. ( |3.7| ) if we use a 
neutrino distribution at zero chemical potential and if we express the field correlators by 
virtue of Eq. ( |3.1U| ) in terms of At,l- 



5 Depolarization in a Relativistic Plasma 

In order to evaluate the depolarization rate explicitly for a specific physical system we 
consider a relativistic QED plasma of the sort encountered in the early universe where the 
chemical potentials of the charged fermions are negligibly small. This system is character- 
ized by the temperature T alone which is taken to be much larger than the electron mass, 
but small enough that muons or pions are essentially absent. 
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Even in such a relatively simple system, the lowest-order polarization functions Ht,l 
depend on K in very complicated ways. In order to arrive at a first estimate we recall that 
Ref. |J implies that the main contribution to the neutrino spin relaxation rate arises from 
spin-flip scattering ui + e ± — > e ± + i/r. The cross section for this process has an infrared 
divergence which is regularized by screening with a scale of order the photon plasma mass 
M. = eT/3. As a first estimate it is thus enough to obtain expressions for Ht,l which 
are precise for ko, k < eT for the leading contribution. In the high-temperature limit 
(m e <C T) they depend only on the variable k /k. Explicitly one finds 



U T (K) 



3M< 



k 2 k 2) 2k [ko-k 




(5.1) 



The real part of these functions is obtained by taking the modulus of the argument of the 
logarithms. Below the light cone (k 2 . < k 2 ) the logarithms yield the imaginary parts 

imn rW = JJf^i-^he^-kl) , 

(5-2) 

Imn L (K) = 3MMi-^ 2 )^o(k 2 -kZ) . 



Above the light cone the only imaginary part of n^L comes from the ie in the retarded 
prescription. Put another way, above the light cone this approximation for Ht,l provides 
support only on the discrete branches which correspond to transverse and longitudinal 
plasmons. 

In the early universe, a typical thermal neutrino momentum is around 3T. Therefore, 
if indeed a typical momentum transfer is of order eT we are in a situation where the 
classical-trajectory approach should be justified where it was assumed that ko, k p and 
also ko T. Therefore, in Eq. ( |4.4j ) the infrared sensitive Bose-Einstein term dominates 
which expands as e(/co)/(e /3 ' fe °' — 1) = T/k . Altogether we thus find for the relaxation rate 
in the classical limit 

r t ot = ^/ dkk dk -—[A T {K)-A L {K)} . (5.3) 

Z7T JO J-k fcn k 



If for the moment we ignore the resummation terms KeIlT,L and Imllr^ in the denominator 
of Eq. (|3.9| ), the classical depolarization rate is found to be 



2rT ,o 2 r°° dk 
i tot = otnT a - I — 
6 Jo k 



(5.4) 
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with a = e 2 /4ir m 1/137 the fine-structure constant. Because this integral diverges it 
needs to be cut off by minimum and maximum momentum transfers, leading to r to t = 
a/x 2 T 3 (2/3) ln(/c max /fc m i n ). Clearly, the infrared cutoff is provided by the plasma mass 
M. = eT/3 while /c max is given by the neutrino momentum itself which is typically 3T. With 
these numbers we find (2/3) ln(/c max //c min ) ~ 2. Because r tot depends only logarithmically 
on the assumed cutoffs we expect this simple estimate to be already rather precise. 

We stress, however, that the initial assumption that only small momentum transfers of 
order eT were important was not justified. In the classical approximation the momentum 
transfers are distributed as 1/k so that the average is (fc) — (fc maj( ^min)/lft(^max/^min)- 
Therefore, with /c max ^> k min the scale for a typical momentum transfer is set by /c max and 
not by fc m j n . Thus, even though the distribution of momentum transfers peaks around eT, 
the average is still of order T and thus not small. 

A numerically precise calculation of r to t thus requires the full quantum expressions. 
Moreover, while the approximate expressions Eqs. ( |5.1| ) and ( |5.2| ) for the polarization func- 
tions are sufficient to implement the screening cutoff at small fco an d k and thus are suf- 
ficiently precise in the denominator of Eq. ( |3.9| ), they are not guaranteed to be accurate 
enough in the numerator of Eq. ( ^.9[ ). While it is well known that the approximate ex- 
pressions are surprisingly accurate even for \K\ of order T, we still use the full one- loop 



expressions which are given in Ref. |15|0. The full one-loop expressions also have the ad- 
vantage of providing continuous support for At,l above the light cone which corresponds 
to pair processes. Therefore, we are able to compare directly the contribution to r to t from 
pair processes (above the light cone) with that from spin-flip transitions (below the light 
cone) . 

The numerical integration of the full quantum expression in Eq. (|4.4|) yields the mo- 
mentum dependent depolarization rate shown in Fig. |I[ When we average these rates over 
a thermal neutrino distribution we obtain 

(r tot ) = 1.8LVT 3 , (5.5) 

where only 1.5% of the coefficient arise from pair processes. They are indeed significantly 
subdominant as predicted in Ref. ||. Moreover, our simple estimate derived from the 
classical approximation was surprisingly accurate. 



6 Early Universe 

6.1 Magnetic Moment Constraint in a Plasma 

The depolarization of the spin in the fluctuating electromagnetic field of the early universe 
affects primordial nucleosynthesis. Although there is a growing awareness that systematic 
errors are large in the determination of the nucleosynthesis limit on the effective number 

3 In Eq. (A. 4) of Ref. Jl5f there is a factor q 2 /u> 2 missing in the term below the light cone. 
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Figure 1: Depolarization rate as a function of neutrino momentum. 



of neutrino species, it still seems reasonable to assume that an extra neutrino degree of 



freedom is not allowed fl6fl . Such an assumption then yields a cosmological limit on the 
neutrino magnetic moment. 

In order to avoid populating the right-handed component of our Dirac neutrinos before 
BBN we need to require that r tot is less than the Hubble rate at all times between the 
muon annihilation and neutrino freeze-out epochs: 



1/2 

, -1 

mpi V 45 



\ ot = 0.0132^ <H = — ^) . (6.1) 



Evidently, the most stringent bound comes from imposing this constraint at as high a 
temperature as possible. We take ~ 10.75 for the effective number of thermal degrees 
of freedom which contribute to the energy density and thus account for the electrons, 3 
left-handed neutrino species, and the photons. With T = 100 MeV this leads to 

fj, < 6.2 x IO^Vb , (6.2) 



where \xb = e/2m e is the Bohr magneton. This result puts previous estimates H, where 
the infrared singularity in the cross section of neutrino-electron scattering mediated by a t- 
channel photon was estimated by a momentum cut-off at the Debye mass, on a solid basis. 
Here we used the full resummed photon propagator to take into account the screening 
effects correctly. 
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To obtain a more precise limit we should examine the relevant Boltzmann equations. 
This is however not really warranted since the most stringent limit on neutrino magnetic 
moments is still obtained from energy loss considerations of helium burning globular cluster 
stars. Plasmon decay would cool these stars too fast unless /i<3x 10~ 12 fi B fTJ\, pi]. Of 



course, this limit is valid only for neutrinos with mass less than a few keV. There are 
additional bounds from SN 1987A [Tj| which are valid up to the experimental v T mass 
limit of about 24 MeV and which are also more restrictive than the cosmo logical limit, 
even though their exact numerical values differ significantly between the works of different 
authors and are based on rather sparse data. Still, unless these astrophysical limits are 
plagued with implausibly large systematic errors we arrive at the conclusion that neutrino 
dipole moments can affect nucleosynthesis only in connection with a large-scale primordial 
magnetic field. 

6.2 Large-scale Magnetic Field 

Another source for spin depolarization would be neutrino interactions with a background 
magnetic field. Let us assume for a moment that such a hypothetical (constant) field exists 
and compare the neutrino spin-flip rate due to a background field to spin flip rate induced 
by the fluctuations of thermal photons. This mechanism has been considered previously 
by several authors |7|, [8], [B| and it is of interest to compare it with the depolarization in 



a stochastic electromagnetic field. In the classical picture the depolarization in a random 
field is a consequence of the random walk that the polarization vector P performs on the 
unit sphere, while in a constant background field the polarization is attenuated by the 
helicity-measuring scattering that projects the coherently rotating polarization vector onto 
the third axis, corresponding to a damping of the off-diagonal elements in the density 
matrix. 

Though it has been derived several times before it easy enough to extract this damping 
directly from Eq. ( |2.3| ). We can simply calculate the eigenvalues of M and see that the 
damping of the third component, in the limit of small D, is 

Au 2 B 2 { D , , 

(6.3) 



At temperatures 1 MeV < T < 100 MeV one finds for effectively massless neutrinos [20] 

that Co>refr = £0°o) with 



e = ^f(. + ^| , (6 .4) 



where C v — \ for electron neutrinos so that £ ~ l.lxlO -20 (T/MeV) 4 , whereas CV = 
for muon and tau neutrinos. The main contribution to D comes from neutrino elastic 
and inelastic scattering with leptons and equals half the total collision frequency |2"T| , p2 . 
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Including scattering with electrons and neutrinos as well as annihilation one finds at tem- 
perature around T = 1 MeV 

D = 2.0AGlT 5 , (6.5) 



for electron neutrinos. Since from Eqs. ( |6.4| ) and ( |6.5D we have that £,(po)/D ~ 100 the 
small D approximation in Eq. ( |6.3|) is good and would remain so even if we include damping 
due to collisions with muons in D. 



Because the electrical conductivity of the universe is large p3[ , a background magnetic 
field is imprinted on the plasma and comoves with the expansion of the universe. Thus 
flux conservation implies that the mean rms field scales with temperature like ^(B 2 ^) 1 / 2 = 
(S 2 ) 1 / 2 = B (T/T ) 2 . We should however point out that in the early universe a large con- 
ductivity implies a large Reynolds number, and hence there is a possibility for turbulence 



which can redistribute magnetic energy to various length scales. This was verified in p4 



where full magnetohydrodynamics was simulated by using a simple MHD generalization of 
the cascade model much used in studying hydrodynamical turbulence. In the magnetic case 
it features a transport of magnetic energy from small length scales to large length scales. 
In a realistic situation the issue is then: has the coherence length of the background field 
grown large enough so it can be treated as a constant mean field on a given scale? Here we 
shall just assume that such a mean field exists. The limit from a large scale background 
field concerns the combination fiB . Requiring that Tf, < H and using Eq. ( |6.3| ) we find 
that 

rp2 /Ox 1/2 f TT vl/2 



T 2 \8J ^^'\D-H 

This constraint should be imposed at the temperature where the right hand side is min- 
imized (and D > H). With the damping rate in Eq. (|6.5| ) and (po) = 3T it happens at 
T = 1.5 MeV which is above the kinetic freeze-out temperature (2D > H gives T > 1 
MeV), and we obtain 

<2.3xl0- 19 ^^ . (6.7) 
MeV 2 V ; 




Here we adopted g* = 10.75 in the Hubble rate. We notice that this mechanism of depolar- 
ization in a large scale field is most efficient at a rather low temperature, in contradistinction 
to the case of a random field which gave a stricter constraint at high temperature. The 
rate T b has been computed previously 0, || but was not applied correctly at the neutrino 
freeze-out in M and an approximate formula for r& was used in |7[] leading to a somewhat 
too stringent bound. Thus, the bound from a large scale field is better than the one from 
small scale fluctuations if B > 3.7xl0~ 9 (MeV) 2 at T = 1 MeV, which is a rather weak 
field compared to other typical scales. The bound on a large scale magnetic field from 
BBN is B < 4 x 10" 5 MeV 2 at T = 0.01 MeV. (If the field is not homogeneous, then 
the limit is less stringent by an order of magnitude.) There is thus room for the large scale 
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field to have been large enough to be the dominating depolarization mechanism without 
contradicting BBN. The existence or non-existence of such a primordial field should be- 
come less speculative with the forthcoming experiments trying to measure the strength of 
the intergalactic magnetic field |26| . 



7 Summary 

We have performed a detailed comparison of the classical and quantum descriptions of 
neutrino depolarization by magnetic moment interactions in the presence of a stochas- 
tic electromagnetic field. Our main result is a general expression of the depolarization 
rate in terms of the electric and magnetic field correlation functions. By virtue of the 
fluctuation-dissipation theorem they are essentially equivalent to the (imaginary parts) of 
the medium's dielectric response functions. Our analysis is exact up to second order in the 
magnetic moment \x and to the order that the electromagnetic field correlation functions 
are calculated. 

We have evaluated the correlation functions and depolarization rate explicitly for the 
case of a relativistic QED plasma. The neutrino depolarization rate is dominated by spin- 
flip scattering on relativistic electrons and positrons. Even though the cross section of this 
process peaks in the forward direction so that the distribution of momentum transfers favors 
values of order eT, an average momentum transfer is of order the neutrino momentum. 
Therefore, the depolarization rate is not approximated well by the classical description 
which is based on the assumption that a typical momentum transfer is small relative to 
the neutrino momentum. A numerically reliable result requires our full quantum treatment. 

Our precise calculation of the depolarization rate puts a previous estimate on a 
firm basis. Imposing the constraint that a right-handed neutrino must not have been in 
equilibrium at nucleosynthesis we derive the limit /i < 6.2 xlO -11 yU# on the neutrino Dirac 
magnetic moment. Other astrophysical limits are more restrictive by about an order of 
magnitude, revealing that neutrino magnetic moments can affect big-bang nucleosynthesis 
only in connection with large-scale primordial magnetic fields which would provide for an 
additional mechanism for left-right transitions. 
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